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Two semi-infinite bodies made of prestressed, homogeneous, Bell- 
constrained, hyperelastic materials are perfectly bonded along a plane 
'"O | interface. The half-spaces have been subjected to finite pure homoge- 

neous predeformations, with distinct stretch ratios but common prin- 
q ■ cipal axes, and such that the interface is a common principal plane of 

strain. Constant loads are applied at infinity to maintain the defor- 
mations and the influence of these loads on the propagation of small- 
amplitude interface (Stoneley) waves is examined. In particular, the 
secular equation is found and necessary and sufficient conditions to 
CN ! be satisfied by the stretch ratios to ensure the existence of such waves 

^& are given. As the loads vary, the Stoneley wave speed varies accord- 

ingly: the upper bound is the 'limiting speed' (given explicitly), be- 
yond which the wave amplitude cannot decay away from the interface; 
the lower bound is zero, where the interface might become unstable. 
The treatment parallels the one followed for the incompressible case 
and the differences due to the Bell constraint are highlighted. Finally, 
the analysis is specialized to specific strain energy densities and to 
the case where the bimaterial is uniformly deformed (that is when the 
stretch ratios for the upper half-space are equal to those for the lower 
half-space.) Numerical results are given for 'simple hyperelastic Bell' 
materials and for 'Bell's empirical model' materials, and compared to 
the results for neo-Hookean incompressible materials. 



1 INTRODUCTION 

This paper is a prolongation of a series of articles concerned with the prop- 
agation of small-amplitude waves in deformed Bell materials. Beatty and 
Hayes [1] wrote the incremental equations of motion for a hyperelastic body 
maintained in a state of static finite homogeneous deformation and subject 
to the constraint of Bell [2J: tr V = 3, where V is the left stretch tensor 
associated with the finite deformation. They studied plane [1] as well as 
torsional |3j waves. Then this author |3] found the secular equation for sur- 
face (Rayleigh) waves propagating on a semi-infinite deformed body made of 
Bell-constrained material. Next [5] the surface stability of such a half-space 
under compression was analyzed and compared with that of incompressible 
rubber. Now the propagation of interface (Stoneley) waves is considered for 
two rigidly bonded semi-infinite bodies made of distinct Bell materials. Each 
half-space was subjected to a finite homogeneous predeformation and the 
principal axes of each deformation are assumed to coincide at the interface, 
which is a common principal plane. Moreover, the interfacial wave is assumed 
to propagate in a principal direction. These common [6j [7] restrictions aside, 
a general necessary and sufficient condition of existence of a Stoneley wave 
is established, together with the corresponding secular equation. 

These general results are obtained in Section 3, after the basic equations of 
the problem have been written in Section 2. For purposes of comparison and 
contrast, these latter equations are written in a manner which is similar to 
those governing the propagation of interfacial waves in deformed hyperelastic 
materials subject to the constraint of incompressibility, det V = 1. Then in 
Section 4, the analysis is specialized to the case where the infinite body 
made of the two bonded half-spaces is uniformly prestrained, that is when 
the principal stretch ratios of strain are the same for the upper half-space 
and for the lower half-space. For "simple hyperelastic Bell" materials, the 
influence of the initial deformation on the wave speed is highlighted. As the 
half-spaces are more and more compressed, this speed tends to zero and the 
secular equation gives the "neutral equation". The critical stretch ratios, 
at which the neutral equation is reached, are obtained. In the case of a 
plane prestrain, these critical stretches, and the corresponding ones for "Bell's 
empirical model" , are compared with the critical stretches obtained by Biot 
[8] for the interfacial stability of (incompressible) neo-Hookean materials, 
often used to model rubber. The conclusion is that, as far as interfacial 
stability is concerned, rubber may be compressed more than Bell's empirical 
materials, but not as much as Bell simple hyperelastic materials, before the 
critical stretches are reached. 



2 PRELIMINARIES 

2.1 Deformed Bell half-spaces 

Consider two semi-infinite bodies made of different homogeneous Bell-con- 
strained hyperelastic materials, separated by plane interface x 2 = 0. Initially, 
each body is at rest in a reference configuration, denoted by Bq for the lower 
half-space x 2 > 0, and by Bq for the upper half-space x 2 < 0. The constitu- 
tive equation for the lower half-space is [2] 

T = j9V + w l + w 2 V 2 , (1) 

expressing the Cauchy stress tensor T in terms of the left stretch tensor 
V. Here p is an undetermined scalar, to be found from the equations of 
motion or of equilibrium, and from the boundary conditions. The constitutive 
parameters ujq and u 2 are defined in terms of the derivatives of the strain 
energy density £ = E(i 2 ,z 3 ) of the hyperelastic body with respect to the 
principal invariants z 2 = [(tr V) 2 — tr (V 2 )]/2, i 3 = det V, as 

u = dZ/di 3 < 0, co 2 = -i 3 ~ 1 dE/di 2 > 0. (2) 

The inequalities above are called the 'Beatty-Hayes ^-inequalities' [2]. Also, 
the Bell constraint is satisfied at all times: 

i x = tr V = 3. (3) 

Similar (starred) equalities and inequalities apply for the upper half-space 
x 2 <0. 

Next, consider that each half-space is subject to a finite static pure ho- 
mogeneous static deformation £> — >■ B\ and Bq — > B* by the application 
of suitable loads. The principal axes of these deformations along x 2 , and 
X\ and 23, two orthogonal directions in the x 2 = plane. Hence the finite 
deformations are described in the coordinate system of the principal axes by 

xr = \rX r (x 2 > 0), and x r = X^X r (x 2 < 0), (T — 1, 2, 3; no sum), 

(4) 
where the Ar, Ap are the principal stretch ratios for each half-space. For the 
lower half-space we have 

V = Diag(A 1 ,A 2 ,A 3 ), 

Ai + A 2 + A 3 = 3, i 2 = A1A2 + A2A3 + A3A1, i 3 = AiA 2 A 3 , (5) 

and the following constant Cauchy stress tensor 

T rr = p X r + u + u 2 \l (r = 1,2, 3; no sum), T tj = (i ^ j), (6) 



where p is constant, and uj , uj 2 are given by (J2J) and valuated at i 2 , 13 given 
by ©3,4, clearly satisfies the equilibrium equations T t jj = 0. So does T*, the 
equivalent starred version for the upper half-space. Finally, the half-spaces 
are rigidly bounded. At the interface x 2 = 0, the continuity of the tractions 
imposes that 

PoX 2 +UJ + X\u 2 =P* \* 2 + UJq +Xfu>2- (7) 

The half-spaces are maintained in these deformed states B±, B\ by the ap- 
plication at infinity of the loads Pr = —T rr , P£ = —T^ r . It follows from (j7|) 
that P 2 = P 2 * • 

2.2 Small-amplitude interface wave 

Now consider that a small amplitude wave travels with speed v and a wave 
number k at the interface, in the direction of x\, with attenuation away from 
x 2 = 0. The following expressions describe this incremental motion B\ — > Bi, 
B\ ^B\: 

x = x + e$l{U(kx 2 )e ik(xi - vt) }, and x = x + e^{V*(kx 2 )e ik{xi ~ vt) }, (8) 

for x 2 > and for x 2 < 0, respectively. Here e is small enough to allow 
linearization and the amplitudes U and U* are of the form U = [Lq, U 2 , 0] T , 
U* = [U*, U2, 0] T . The remainder of this paper is devoted to the study of 
these interfacial (Stoneley) waves. Dowaikh and Ogden [6] and Chadwick [7] 
conducted similar studies for prestrained incompressible and compressible 
materials, respectively. In both cases, the authors proved that the initial 
normal traction T 2 2 = T 2 * 2 played no role in the resolution and analysis of 
the problem. It can be checked that such is also the case for deformed Bell 
materials and without loss of generality, we choose T 2 2 = T 2 * 2 = 0, and thus 
P 2 = P 2 * = (of course, the traction T22 might play a role in issues other 
than the derivation of the secular equation for Stoneley waves.) Hence the 
following loads Pi = — T n , P 3 = — T 33 and P*, P 3 * are applied at x± = ±00 
and x 3 = ±00 in order to sustain the finite homogeneous deformation, 

Pr = (A 2 - A r )(-o;o + ArAsc^)/^, 
Pr = (A2 ~~ A r )(— Uq + X r X 2 u 2 )/X 2 , 

(r = i,3). (9) 

Moreover, this assumption enables us to use previously established results for 
the incremental equations of motion. Specifically, the equations are written 
as the following system of first order differential equations for the components 



of the displacements U, U*, and of the tractions t, t* on the planes x 2 = 
const, (see Ref.JI] for details): 

b 3 \ 2 2 U[ + ib 3 \ 2 2 U 2 - £1 = 0, t[ + zAiAj 1 ^ - (AiAj^C - pv 2 )XJ x = 0, 
U' 2 + zAxA^ 1 ^! = 0, t' 2 + ih - [b 3 (Xl - Xj) - pv 2 ]U 2 = 0, (10) 

where 

— UJq + Ai\ 2 LiJ 2 

3= A 2 (A 1 + A 2 ) ' 

C = A 1 A2C11 + AiA 2 C22 — C\2 ~~ C21 — 2o;o — (A 1 + X 2 )uj 2 , 
Cij = 2X 2 i S ij uj 2 - Aj(w 02 + A-w 2 2) + AiA 2 A 3 (u;o3 + A?w 23 ), (11) 

for x 2 > 0, and a starred version for x 2 < 0. Here p (p*) is the mass density 
of the lower (upper) half-space, and the derivatives u r, <^2r (r = 2, 3) of the 
material parameters ojq, uj 2 are taken with respect to ir and evaluated at i 2 , 
is given by (jSJ). The half-spaces are rigidly bonded at x 2 = so that the 
following boundary conditions apply, 

U(0) = U*(0), t(0)=t*(0), (12) 

together with the requirement that the wave vanishes as x 2 — > ±00. 

3 GENERAL PREDEFORMATIONS AND 
STRAIN ENERGY FUNCTIONS 

Here we solve the problem at hand, that is we solve the equations of motion 
(HOP and then apply the boundary conditions (j!2p in order to derive the 
secular equation for the speed of Stoneley waves in deformed Bell materials. 
We adopt an approach reminiscent of that used by Dowaikh and Ogden [6] 
for Stoneley waves in incompressible materials and show how each constraint 
leads to different results. 

3.1 Secular equation; comparison with incompressible 
materials 

For small-amplitude waves of the form (jHJ)i, the incremental constraint of 
incompressibility, u^i + u 2i2 + 11,^3 = 0, imposes that 

iU x + U 2 = 0, (13) 



suggesting the introduction of the function (p defined by 

U x {kx 2 ) — i(p'{kx 2 ), U 2 (kx 2 ) — <p(kx 2 ). (14) 

In our context, the incremental constraint of Bell, Aiii^i+A^^ 2+^3^3,3 = 
0, imposes Eq. ffTUjU . that is 

iXxX^Ux + U 2 = 0, 

suggesting the function tp defined by 

Ui{kx 2 )=iy\k\i\ 2 ~ 1 x 2 ), U 2 {kx 2 ) = <p>{k\ 1 X 2 l x 2 ). (15) 

From (TXOTl xi 1. the traction components are also expressed in terms of (p and 
its derivatives: 

h = ib z Xl(X 1 X^"+^), t 2 = -b 3 X 1 X 2V "'+X^ 1 X 2 (X 1 X 2 l C-hX 1 X 2 -pv 2 ) v '. 

(16) 
Finally, Eg. (TIT)]) 9 reads 

b z X\<p"" - (X x X 2 l C - 26 3 AiA 2 - pv 2 )ip" + (b 3 XJ - pv 2 )ip = 0. (17) 

The same procedure is of course also valid for the upper half-space, with the 
introduction of a function ip*. 

Dowaikh and Ogden [6] expressed the strain energy density as a function 
W(Xi, X 2 , A3) of the principal stretches of deformation, rather than as a func- 
tion of the invariants E(ii,i 2 , 13), a practice also favored by Biot [9]. With 
this choice, and the introduction of the function (p in ffT5|) . they proved that 
for incompressible materials, the counterpart to (fT7|) is 

7^"" - (2/3 - pv 2 )(p" + (&- pv 2 )<p = 0, (18) 

where a, ft, and 7 are defined in terms of W and its derivatives with respect 
to the Aj as 

aX 2 2 = 7 A? = (A1W1 - X 2 W 2 )XiXj/(Xl - X 2 2 ), 

2(3 + 2 7 = XJWn + X\W 22 - 2X X X 2 W 12 + 2X 2 W 2 . (19) 

Also, for incompressible materials, the assumption of strong ellipticity for 
the equations of motion implies that 

d>0, 7 > 0, (3+^/d^>0. (20) 



For Bell constrained materials, introduce by analogy a, (3, and 7, defined 
by 

a = b 3 Xl 7 = Ml, 2/3 + 2^ = AiA^C, (21) 

or equivalently, by 

a\ 2 2 = 7 A 2 = (Wx - ^ 2 )A 2 A 2 /[A 3 (A 2 - X 2 2 )], 

2/3 + 2^ = (W11 + iy 22 - 2Wi 2 )Ai/(A 2 A 3 ). (22) 



Now rewrite the equation of motion ( 1171) as 

V" - (2/3 - pv 2 )tp" + (« - P^ 2 V = 0. (23) 

Using the strong ellipticity condition for the equations of motion in deformed 
Bell materials (see Appendix) together with Eqs.([2]) and ( TTTj) . we also find 
the following inequalities, 

a > 0, 7 > 0, /3 + a > 0. (24) 

The contrast between the effects of incompressibility and of the Bell con- 
straint is striking (compare for instance the first term in (TT8|) and in (|23j) .) 
Now Equation (|2"31 and its starred version for the upper half-space are solved. 
Because Stoneley waves vanish with distance from the interface, the so- 
lutions (p and ip* must be of the form 

<p(z) = Ae- SlZ + Be- S2Z {z > 0), <p*(z) = A* e s i z + B* e s * z {z < 0), (25) 

where si 7^ s 2 , 3ft(sj) > 0, s* 7^ s*>, 9ft(s.*) > 0, for some constants A, B, A*, 
B*. Explicitly, the Sj, s* are roots of the biquadratics 

as 4 - (2/3 - pv 2 )s 2 + a - pv 2 = 0, aV 4 - (2/3* - pV)s* 2 + a* - p*v 2 = 0, 
s 2 + s 2 = (2/3 - pv 2 )/a, sf + sf = (2/3* - p*v 2 )/a*, 

2 2/ 2\ / *2 *2 / * * 2\ / * 

5^2 = (a — pv J/a, 5 1 5 2 = ( a "~ P v J/ a • 

(26) 

The roots sf and s 2 . of the real quadratic (l2"6"lb are either both real or both 
complex; if they are real, they are non-negative, otherwise si and s 2 are 
purely imaginary in contradiction with the decaying requirement 9ft(sj) > 0; 
if they are complex, they are conjugate, because the coefficients of (12"B"1) i are 
real. In both cases, s\s\ > and similarly, s* 2 s*, 2 > 0, so that 



< v < min{t>L,f2}, where Vl = y/a/p, v* L = \Ja* / p* . (27) 



In fact, the interval of possible values for v for which the wave decays away 
from the interface (the subsonic interval) might be smaller than the interval 
defined above. This point is clarified in §3.21 For the time being, we derive 
an explicit expression for the secular equation. Equations (TT21 express the 
continuity of the displacement and traction incremental amplitudes at the 
interface x 2 = 0. Using the expressions ([Hj) , (TTB|) of C/j, U*, t i: t*, in terms 
of tp, tp* and their derivatives, they lead to 

(p(0) = <p*(0), ^(0) = tp*'(0), 

7k/V(0) + ?(0)] = 1*[\^*"(0) + <P*(0)}, (28) 

V a \ a* 



a 7 y/"(0) + \l±(2p + y/arj - P« V(0) = 

V^tV'W + \ — (2/3* + v^T - P^ 2 )</(0). 



f ,-y* 

a* 

Using (J2B]) q/i and (j2"5j) . we rewrite this system as four homogeneous linear 
equations for the unknowns A, B, A*, B*\ 

A + B = A* + B\ 

s ± A + s 2 B = -s\A* - s* 2 B*, 

(^forjsl + 7) A + {^/oqs\ + l)B = 

(y/aFfs? + -f*)A* + (y^Ffsf + 7*)fi*, 
Si{y/oPys\ + 7)^ + s 2 (a/o7Si + 7)5 = 



s\{y/cFfsf + 7*)A* - s^a/o^Ts;; 2 + 1*)B*. (29) 



The vanishing of the corresponding determinant yields the secular equation. 
Dropping the factor (si — s 2 )(s^ — s 2 ) and using the quantities 77, r (and 77*, 
r*) defined as follows, 



n=\ ?L -^, r=J{T, + l)* + 2?—Z t (30) 

a V a 



the secular equation is written compactly as 

• / \ 2r a 2 / a 1 \2l 1 *2r Q: * *2 / / a * i\2i 



^ r /aa* . ^ N . /a . . /a* 

+ 77* \/ — fa + ^ rr + 2 \/-^ " X V T ^ " X 31 

V 77 V 7 V 7 

= 0. 

The quantities rj, rf are real by fl27|) and so are r, r*, at least as long as the 
decaying condition is satisfied. We now elaborate on this last point. 
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3.2 The subsonic interval 

From ( 126]) and the definition (130|) of r, expressions for the squared sum and 
difference of the roots Si follow: 

(s 1 + s 2 ) 2 = r 2 , ( Sl -s 2 ) 2 = s 2 , (32) 

where 



■2 2/3 — pv at — pv 2 2 8 — a ,, 



s 2 = ' r - 2\ ^—^- = ( v - l) 2 + 2'—^ = r 2 - 477. (33) 

a V ex a 

Hence the requirement of exponential decay 9R(sj) > is equivalent to 

9fc(r ± s) > 0. (34) 

Now we span the subsonic range, that is the interval / (say) of possible values 
for v such that (1341) is satisfied. This interval is I = [0, v] where v is called 
the limiting speed. 

Starting at v = 0, we have, using the strong ellipticity condition (|2"4"]U . 



v = n = 1, r---- \l2'—^ > 0, r 2 - s 2 = 4 > 0, (35) 

a 



so that (1541) is satisfied. 



Now, if 2/3 — a > 0, then we may increase v from to v l = y ac/p, where 



v = VL n = 0, r == W — — - > 0, r 2 - s 2 = 0, (36) 

V a 

and ?7, r, r 2 — s 2 , decrease monotonically with v but remain non-negative 
real numbers, so that (1341) is satisfied over the interval [0, Vl\. Clearly in this 
situation, the limiting speed is v = v l = yotjp. 

However, in the situation where 2/3 — a < 0, we may increase v from 
only up to v (where r = 0) in order to satisfy (134")) . where v is defined by 



pv 2 = 2[/3 - a + y / 2a{a - 3)} < a = pv\. (37) 

Clearly the limiting speed is then v = v. 

Conducting a similar analysis for the upper half-space, we conclude that 
the limiting speed is 

mm{v L ,v* L }, when 2/3 > a, 23* > a*, 

min{{;,t;2}, when 2/3 < a, 2/3* > a*, , . 

mm{v L ,v*}, when 2/3 > a, 2/3* < a*, ^ ' 

min{-u,{;*}, when 2/3 < a, 2/3* < a*. 
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Note that the analysis conducted in this subsection (defining the subsonic 
interval) and in the following one (defining the conditions of existence of a 
Stoneley wave) rely directly upon the methods developed by Chadwick [7j 
for compressible materials, with the modifications required to accommodate 
the Bell constraint. 



3.3 Existence and uniqueness of interfacial (Stoneley) 
waves; comparison with surface (Rayleigh) waves 

Here we derive the conditions of existence for a Stoneley wave at the interface 
of two deformed Bell-constrained half-spaces, using a "matrix reformulation" 
of the secular equation (l3Ti) . a method suggested by Chadwick [7] and based 
on the surface impedance method of Barnett et al. [10]. Indeed, introducing 
the following symmetric 2x2 matrices 



M(v) 



7 



rjr 1 — 



n 



M*(v) 



7 



7* ' 

'-v* - 1 



%ri* - 1 

7* ' 



7* 



we find that the secular equation (131(1 corresponds to 



i(v) 



det N(v) = 0, where N(v) = M(v) + M*(v) 



(39) 



(40) 



From then on, it is an easy matter to transpose Chadwick's results [7j for 
incompressible materials to Bell materials, and to show inter alia that the 
eigenvalues of N(v) decrease monotonically as v increases in /. Then the 
following results apply (see also Barnett et al. [10] for Stoneley waves in 
linear anisotropic elasticity): when an interfacial wave exists, it is unique; it 
propagates at a speed which is greater than the speed of the Raleigh wave 
associated with either each half-space; it exists if and only if 



i(0) > 0, i{v) < 0. 



(41) 



In the (Ai, A2, A3, A*, A2, A^-space, the curve i(v) = is called the limiting 
equation and the curve i(0) = is the neutral equation. Dowaikh and Og- 
den [S] refer to this latter equation as the "exclusion equation" because it 
rules out the possibility of incremental inhomogeneous static deformations 
in homogeneously deformed half-spaces. Biot [8] called it the "characteristic 
equation for instability" , because it "corresponds to the spontaneous appear- 
ance of sinusoidal deformations at the interface." We obtain this equation 
for Bell materials by taking 77 = rf — 1 in the secular equation (|3T|) : 



V / ^(M + V / 2rP^±(^-v / aY + f-7) = 0. (42) 
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4 UNIFORM PREDEFORMATION AND 

SPECIFIC STRAIN ENERGY FUNCTIONS 

In the preceding section, we covered some ground on the propagation of 
Stoneley waves in deformed Bell materials. We obtained results for general 
prestrains and unrestricted strain energy density functions. We now turn our 
attention to special configurations. First, we assume that the infinite body 
made of the two bonded Bell-constrained semi-infinite bodies is predeformed 
uniformly in the whole space, so that the principal stretches are the same for 
the lower and upper half-spaces: 

AJ = Ai, A^ = A 2 , A* = A 3 . (43) 

This deformation is possible with the application of the following loads, 

^r = (A 2 - A r )(-Wo + A r A 2 w 2 )/A 2 , 
p* = (A 2 - X r )(-tu* + A r A 2 w 2 *)/A 2 , 
(r = l,3). (44) 

Next, we specialize the analysis to specific forms of the strain energy 
functions and make the connection with historical results obtained for in- 
compressible materials. 

4.1 Simple hyperelastic Bell materials 

Here we consider that the lower and upper half-spaces are made of "Simple 
hyperelastic Bell" materials, for which the strain energy functions are of the 
form [2], 

Wshb = C\(3 — AiA 2 — A 2 A 3 — A3A1) + C 2 (l — AiA 2 A 3 ), 

Wshb = CJ(3 - A X A 2 - A 2 A 3 - A 3 A 3 ) + C*(l - AxA 2 A 3 ). (45) 

The material constants Cj, C* satisfy, according to the ^-inequalities (J2J), 

d > 0, C 2 > 0, C{ > 0, C* 2 > 0. (46) 

Now, the expressions (I22|) and (J5D] )? for a, (3, 7, and r (and for a*, (3*, 7*, 
and r*) reduce to 

7 = (dA^ 1 + C 2 )A 2 /(A 1 + A 2 ) > 0, a = (3 = A?AJ 2 7 > 0, r = v +l. (47) 
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Because here 2/3 — a = a > 0, we deduce from the discussion undertaken in 
§3. 21 that the limiting speed for two bonded simple hyperelastic Bell materials 
is 

v = min{ya7/7, a/«*/7*}. (48) 

For v G I = [0,v], the secular equation (|3T|) reduces to 

\^ 2 \ 2 Av) ee 7 V + V 2 + (1 + 2Ar 1 A 2 )r? - A^A*] 

+ 7*V 3 + ^/* 2 + (1 + 2Ar 1 A 2 )?f - Ar 2 A"] (49) 

+ TfKv + V*)(V + 1W + 1) + 2(Ar : A 2 - r/)(Ar 1 A 2 - r/*)] 
= 0. 

The conditions of existence of a Stoneley wave are (14T1) . with the above 
simplifications. Explicitly, the neutral condition i(0) > factorizes to 

[(3 7 * + 7)Ai - (7* - 7)A 2 ] [(37 + 7 *)Ai - (7 - 7*)A 2 ] > 0. (50) 

The sign of each factor depends on the sign of 7 — 7*. Introducing the 
quantity e defined by 



= (Ci+C 2 A3)/(Q + C 2 *A 3 ) when 7 < 7*, 
7 77=(C 1 * + QA 3 )/(C 1 +C 2 A 3 ) when 7* < 7, l ° iJ 



we write the neutral condition (|50l) as 

(3 + e)Ai - (1 - e)A 2 > 0. (52) 

Note that the neutral condition (1521) gives an implicit relationship between 
the stretch ratios Ai and A 2 because A3 appearing in e is linked to Ai and 
A 2 through the Bell constraint: A 3 = 3 — Ai — A 2 (explicitly, the relationship 
between Ai and A 2 is quadratic.) However, for the plane strain A 3 = 1 (and 
A 2 = 2 — Ai), e is independent of the stretch ratios and the critical stretch 
(Ai) cr , at which the neutral equation is reached, is 

(Ai) cr = ±(l-e). (53) 

Note also that when one half-space is absent (7 = or 7* = 0) then e = 
in (152]) and we recover the relative universal surface stability condition 
3Ai — A 2 > for the remaining half-space [5]. 

Finally we express the limiting condition i(v) > in the case where 
v l = \fot*7~P* < v l — yoijp- Then we have 




a */-\ * r, t*s a-pvj 1 pa" 



,*2 



p* V a v P a 

(54) 
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and the limiting condition i(v) > is 

Vl(vl + l)( V L + l + 777) - fa ~ Ar 1 A 2 (l - 7*/7)] 2 < 0. (55) 

In the case where Vi = sfotfp < v* L = ^Ja* / ' p*, the starred and unstarred 
quantities above are interverted. 

In his seminal paper, Stoneley [IT] showed, by means of numerical exam- 
ples, that the existence of 'waves at the surface of separation of two solids' 
depended heavily on the material properties of each half-space. We treat 
two numerical examples below, choosing p, a, 7, p*, a*, 7*, such that a 
connection is made with his results. 

Example 1: a*/p* = a/p, with p = 8.2, p* = 3.2. 
In this case, we have 

7* a* p* 3.2 . , 

e = — = — = — = — = 0.39024. 56 

7 a p 8.2 v ; 

Important simplifications occur in this special case, as 



pv 2 / p*v 2 



a V of 

a 




v = v L = v* L 

" p 

Vl = V*l = 0. (57) 

Hence, the limiting condition i(v) > is automatically satisfied, because 
T)l = in (|55p . The only condition of existence of a Stoneley wave is the 
neutral condition ([52]), that is 



Ai-0.17986A 2 >0. (58) 

Finally, the secular equation (J4"9l reduces further to 

(1 + e)V + V 2 + (1 + 2Ar 1 A 2 )7? - Af 2 ^] + 4e( V - \^\ 2 ) 2 = 0, (59) 

with e given in (1561) . In particular, when the materials are undeformed in 
the static state (Ai = A 2 = A3 = 1), the constraint of Bell coincides with 
the constraint of incompressibility [2] , and we recover Stoneley's result [TTJ : 
v = 0.99287£, for linear isotropic incompressible materials with material 
parameters matching those of this example. 

Figure 1(a) shows the intersection of the triangle Ai + A 2 + A3 = 3 of pos- 
sible values for the stretch ratios in Bell materials with the neutral condition 
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(a) Region of stability 



(b) Stoneley wave speed 



Figure 1: Simple hyperelastic Bell materials (a*/p* = oc/p, p = 8.2, p* 
3.2). 



(158"|) for the combination of simple hyperelastic Bell materials in this example; 
the visible part of the triangle represents the configurations where a Stoneley 
wave exists, and where the interface is stable with respect to incremental 
deformations. 

Figure 1(b) illustrates the influence of the loads on the Stoneley wave 
speed (thick curve) in the case of plane strain A 3 = 1. Then A 2 = 2 — Ai, 
and the critical load obtained from flo'81 is (Ai) cr = 0.30488. The speed, 
scaled with respect to the limiting wave speed v = \fa/p = \/ot* j ' p*, is the 
coordinate on the vertical axis; the stretch ratio Ai is the coordinate on the 
horizontal axis. In tension (1 < Ai < 2) the Stoneley wave propagates at a 
speed which is within less than 1% of the limiting speed; under increasing 
compression (1 > Ai > (Ai) cr ), the speed drops rapidly to zero. The Figure 
also shows the scaled Rayleigh wave speed associated with either half-space 
[5] , computed by taking e = in (1591) (thin curve, always below the Stoneley 
wave speed curve as expected). 

Example 2: a*/p* = \aj ' p, with p = 2, p* = 1. 



In this case, we have 



The limiting speed is 



r 



a* 3p* 



7 a 2p 4 



v = v L 




(60) 



(61) 
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and the corresponding 77's are 










Vl = 0, rf L = y 1 - 


a* 


/- 


p*a 
pa* 


1 



(62) 

The limiting condition (155]) (with the starred and unstarred quantities 
interverted) is now satisfied for Ai, A 2 such that 



(^10 + 8^- v / 3)Ai-A 2 <0, or Ai - 0.31723A 2 < 0. (63) 

Here we see that no interface wave may propagate when the materials are 
undeformed in the static state (Ai = A 2 = A3 = 1), as noted by Stoneley 
[TT] for linear isotropic incompressible materials with material parameters 
matching those of this example. 

The other condition of existence of a Stoneley wave is the neutral condi- 
tion (152]) . that is here 

15Ai - A 2 > 0. (64) 

Figure 2(a) shows the intersection of the triangle Ai + A 2 + A 3 = 3 of 
possible values for the stretch ratios in Bell materials with the limiting con- 
dition (1631) and with the neutral condition (I64p for the combination of simple 
hyperelastic Bell materials in this example; the plane of the Figure coincides 
with the plane of the triangle. We see that the range of possible stretch ratios 
is far smaller here than in the previous example. For instance, in the case of 
plane strain A 3 = 1, the stretch ratio Ai must belong to the (compressive) in- 
terval ]0.125, 0.48167[. Note that in this plane strain case, the Stoneley wave 
exists in a range where Rayleigh waves do not exist (in [5] it is proved that 
the domain of existence of Rayleigh waves in plane strain is: Ai GJ0.5, 2[.) 

Figure 2(b) illustrates the influence of the loads on the Stoneley wave 
speed (thick curve) in the case of plane strain A3 = 1. The stretch ratio Ai 
is the coordinate on the horizontal axis. The speed, scaled with respect to 
the limiting wave speed v = \/a/p, is the coordinate on the vertical axis; 
it is computed by solving numerically PUj) for y = v/v with i] = a/1 — y 2 , 

rf = \ 1 — \y 2 -, and 7* = I7. Also represented are the Rayleigh wave 

speeds (thin curves) associated with each half-space, computed by solving 
numerically (T4T?1) for y = v/v with 7] = \Jl — y 2 and 7* = (lower thin 



curve, with y = 1 as an horizontal asymptote) or?/* = wl — |y 2 and 7 = 
(upper thin curve, with y — * /| as an horizontal asymptote). 



15 





(a) Region of stability 



(b) Stoneley wave speed 



Figure 2: Simple hyperelastic Bell materials (a*/p* = ^a/p, p = 2, p* = 1). 

4.2 Bell's empirical model 

For Bell's empirical model materials [2], the strain energy functions of each 
half-space are 



Wbem = |^o[2(3 - AiA 2 + A 2 A 3 + A 3 Ai)]4, 



IF, 



BEM 



z/ *[2(3-A 1 A 2 + A 2 A 3 + A 3 A 1 )]4 : 



2 -, ,* 
3 



(65) 



where z/o, v§ are positive constants. For these models, the material response 
functions ujq, ^2, ^0 and ui 2 provided by ([2]) are 



u = u* Q = 0, u 2 = i 3 %[2(3 - i 2 )] 4, uj* = € u 2 , e = Uq/u q , 
with i 2 , 13, given by ©3,4- Then a, (3, 7, a*, (3*, 7* are computed as: 
Xi A x A x r (Ai — A 2 ) 1 A X A 2 



A1 + A2 

a* = ea, (3 



UJ 2 



X 



2 7, 



£ 



4(3 - z 2 ) 



w 2 - 



A1 + A2 



w 2 , 



e/3, 7* = £7. 



(66) 

(67) 
(68) 



Although the strain energy function of a "Bell's empirical model" material 
depends only upon one material constant (z/ ), as opposed to two (C\ and 
C 2 ) for a "simple hyperelastic Bell" material, the full analysis of the Stoneley 
wave existence is too cumbersome and lengthy to be followed here. Indeed, 
we saw in §3.21 that the definition of the limiting speed depends on the sign 
of 2/3 — a which, for these materials, is equal to 

2/3 -a = oo 2 Xl[3Xl - (12 - 5A 2 )Ai(Ai + A 2 ) - A|]/[4(3 - i 2 )(Xi + A 2 )]. (69) 
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The sign of this quantity (and hence the definition of v) depends heavily 
on the values of the stretches Ai and A2 (For instance in the equibiaxial 
case Ai = A3, we have 2(3 — a = w 2 A 2 (7Ai — 9)/[4(3 — Ai)], which changes 
sign at Ai = |.) Consequently, the limiting condition is difficult to obtain 
in general. Moreover, the connection with Stoneley's results [11 J cannot 
be made because the parameters u 2 , u> 2 are singular when the material is 
isotropic (then i 2 — 3), where a different approach must be adopted |2J. 

Nevertheless, the neutral equation, which is independent of v and takes 
place away from isotropy, can be obtained and compared with Biot's "stabil- 
ity equation" for rubber-like incompressible materials [8]. We deduce it by 
specializing f|42|) to the values ( 1ST]) of a, (3, 7, a*, (3*, 7*, as 

4Ai (Ai-A 2 ) 2 rl-ei 2 (Ai-A 2 ) 2 



Ai + A 2 4(3 - i 2 
In the case of plane strain A3 



Ll 



Ai(AiH-A 2 ) 



0. 



(2Ai - l)Aj 



1, the equation reduces to 
1- 



(Ai 



1 



0. 



(70) 



(71) 



Note that when e = in (I70|l . only the lower half-space subsists and we re- 
cover the surface bifurcation criterion for Bell's empirical materials in com- 
pression [T2] , 

(Ai-A 2 ) 



X 1 A 2 



0. 



(72) 



4(3 -i 2 ) 

Figure 3 (where the plane of the Figure coincides with the plane of the 
triangle) shows the intersection between the triangle of possible stretch ratios 
for Bell materials and the neutral curve in the cases where e = 0.0 [5] (thickest 
curve), e = 0.2, e = 0.4, e = 0.6 (thinnest curve), and e = 0.8 (dotted curve). 



4.3 Comparison with neo-Hookean incompressible ma- 
terials 

Biot [8] investigated the surface instability of incompressible materials under 
finite compression, and obtained the neutral condition, which he called the 
"characteristic equation for interfacial instability." He then discussed the 
equation in the particular case of a neo-Hookean strain energy function, 

W nH = /io(A? + A 2 + A 2 -3), (73) 

and expressed it, in the case of plane strain A3 = 1, as: 

-l + A 2 i 2 rl-^ 2 



K 



Ll- AfJ 



Ll + eJ 



^0 
0o" 



(74) 
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Figure 3: Neutral curves for 'Bell empirical model' materials. 

He did not mention that the result was also valid for the neo-Hookean strain 
energy function, 

WW = V 1 (\j + \j + A^ - 3) + V 2 (\l\l + X 2 2 Xl + \ 2 3 \j - 3), (75) 

by replacing the rigidity ratio e with 

T)* — T)* 

< = %%■ (76) 

Table 1 shows the numerical values of the critical stretches for the classes 
of Bell's empirical model (2nd column), of neo-Hookean and Mooney-Rivlin 
incompressible materials (3rd column), and of simple hyperelastic Bell ma- 
terials (4th column) in the case of plane strain, for different values of the 
rigidity ratio e. It appears that rubber can be compressed more than Bell's 
empirical model but less than simple hyperelastic Bell materials, before the 
neutral equation is satisfied. 

Table 1: Critical stretch ratios for interface instability (A3 = 1) 



e 


Bell empirical 


rubber 


simple Bell 


0.0 


0.6667 


0.5437 


0.5000 


0.2 


0.6105 


0.4457 


0.4000 


0.4 


0.5625 


0.3393 


0.3000 


0.6 


0.5266 


0.2257 


0.2000 


0.8 


0.5060 


0.1085 


0.1000 
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Appendix 

Strong ellipticity conditions for the incremental equations of motion in Bell 
materials. 

For an unconstrained hyperelastic material maintained in a static state 
of pure homogeneous deformation, the incremental equations of motion are 

AjUkUk,ji = puijt, (1) 

where A. is the fourth-order instantaneous linear elasticity tensor. The 
nonzero components of A. are: 

Aijij = (W, - A^)A 2 /(A 2 - A, 2 ), (2) 

(no sum) when the underlying deformation has distinct principal stretch 
ratios Ai, A2, A3. The assumption of strong ellipticity for the equations of 
motion ([!]) imposes that 

•Ajiikvrijnimink > 0, for all nonzero m, n. (3) 

When the material is incompressible, det V = 1, where V is the left 
stretch tensor, and an arbitrary pressure P is introduced [13J. Then the 
equations of motion are strongly elliptic for all values of P when the inequal- 
ities (|H]) hold subject to 

m ■ n = 0. (4) 

When the material is subject to the Bell constraint tr V = 3, an arbitrary 
scalar p is introduced [2] . Then the equations of motion are strongly elliptic 
for all values of p when the inequalities (J3j) hold subject to 



m • Vn = 0. (5) 



The unit vectors 



1 _i _i _i 



m = X t cos d'x + A 2 sin0j, n = — X 1 sin 9i + A 2 cos0j, < 6 < 2ir, 

(6) 
are two such vectors, and the strong ellipticity condition says that 

A cos 4 6 + 2B sin 2 9 cos 2 9 + C sin 4 9 > 0, (7) 
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for all 9, where A, B, C, are given by 

A — X 1 X 2 A\2\2i C — X 1 A 2 ^2121, 

B = (ArMmi + X 2 2 A 2 222)/2 - X^X 2 l (A ll22 + A l221 ). 



Choosing 9 = arctan(A/C) 4 ; we arrive at B + y/AC > 0, which is equivalent 
to 

(3 + a>0, (9) 

where 

a = (W l -W 2 )X 2 1 /[X 2 Xz{X\-Xl)} ) 

(3 = (W u + W 22 - 2W 12 )A 1 /(2A 2 A 3 ) - aA 2 /Ai. (10) 



21 



